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Effect of exchange anisotropy on the relaxation time of spin and vector chirality is stud- 
ied for the antiferromagnetic classical Heisenberg model on the triangular lattice by using 
the nonequilibrium relaxation Monte Carlo method. We identify the Berezinskii-Kosterlitz- 
Thoulcss (BKT) transition and the chiral transition in a wide range of the anisotropy, even for 
very small anisotropy of ~ 0.01%. As the anisotropy decreases, both the critical temperatures 
steeply decrease, while the BKT critical region becomes divergently wide. We elucidate a sharp 
"V shape" of the phase diagram around the isotropic Heisenberg point, which suggests that 
the isotropic case is exceptionally singular and the associated Z2 vortex transition will be iso- 
lated from the BKT and chiral transitions. We discuss the relevance of our results to peculiar 
behavior of the spin relaxation time observed experimentally in triangular antiferromagnets. 

KEYWORDS: geometrical frustration, triangular lattice, Heisenberg model, nonequilibrium relaxation, 
BKT transition, chiral transition, Z2 vortex transition 
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Antiferromagnet on the two-dimensional triangular 
lattice has been intensively studied as one of the most 
fundamental models for the geometrically frustrated sys- 
tems. 1 For the isotropic Heisenberg model with nearest- 
neighbor interactions, it is believed that the ground state 
of the system exhibits a thrcc-sublatticc 120° long-range 
order, 2 whereas the magnetic ordering is no longer re- 
tained against thermal fluctuations. 3 Nevertheless, an 
interesting possibility was proposed by Kawamura and 
Miyashita, 4 ' 5 that is, an unconventional topological tran- 
sition at a finite temperature (T) — Z2 vortex transi- 
tion. From the symmetry point of view, the Z2 vortex 
transition is different from the conventional Bcrczinskii- 
Kosterlitz-Thouless (BKT) transition which occurs in 
the presence of anisotropy. 6. The relation between these 
two topological transitions, however, is not fully under- 
stood yet. In particular, it is still unclear how the system 
behaves in the region of vanishing anisotropy. 

Experimentally, several materials with triangular lay- 
ered structure have been studied, and recently, the Z2 
vortex attracts renewed interests for understanding of 
their peculiar properties. One of the peculiar proper- 
ties is anomalous enhancement of the spin relaxation 
time. Critical divergence of the relaxation time is ob- 
served in an anomalously wide range of T in many com- 
pounds, such as ACv0 2 (^4=Li,H,Na), 8 ~ n Li 7 Ru0 6 , 12 
and NiGa2S4. 13, 14 The critical behaviors are often ar- 
gued to be a fingerprint of the Z2 vortex transition. The 
Z2 vortex, however, is a topological object specific to 
spin-rotational-invariant systems, and hence, it is not 
trivial whether its influence is observed in real com- 
pounds in which anisotropy exists. 

In this letter, to shed light on the origin of the anoma- 
lous critical behavior and its relation to the Z2 vortex 
transition, we directly calculate the relaxation time in the 
antiferromagnetic Heisenberg model with classical spins 
on the triangular lattice. We focus on how the anisotropy 



in exchange interactions affects the behavior of relaxation 
time. We determine the finite- T phase diagram precisely 
by varying the anisotropy, and uncover the exceptionally 
singular nature of the isotropic Heisenberg case. 
Our model Hamiltonian is defined in the form 



\Sj Sj ), 



(1) 



where J is the antiferromagnetic exchange interaction, 
and Si = (Sf , Sf , Sf) is a vector representing the clas- 
sical spin at site i (we normalize as \Si\ = 1); the 
summation (ij) runs over the nearest-neighbor bonds of 
the triangular lattice. We introduce here the exchange 
anisotropy A. For the XY anisotropy (A < 1), it is known 
that both BKT and chiral transitions occur at different 
but very close temperatures. 15 ' 16 On the other hand, for 
the Ising anisotropy (A > 1), it is known that two differ- 
ent BKT transitions occur separately for the longitudi- 
nal S z and the transverse (S x , S v ) components. 17, 18 It is, 
however, still controversial how these four transitions be- 
have as the system approaches the isotropic Heisenberg 
point A = 1. We will discuss this issue later. 

We calculate the relaxation time of the model (1) by 
using the nonequilibrium relaxation (NER) method. 19 In 
this method, the relaxation time is directly computed by 
analyzing the relaxation process from an initial ordered 
state in terms of the Monte Carlo (MC) dynamics. We 
typically perform the relaxation up to 10 5 MC steps by 
using the standard Metropolis local update for the sys- 
tem size N s = L x L up to L = 4002 under the periodic 
boundary conditions. We confirm that the finite-size ef- 
fect is negligibly small. The results are averaged over 
eight independent MC runs. We choose the initial state 
to be a three-sublattice 120° state; spins are set to be in 
the xy plane for A < 1, while they are in the xz plane with 
aligning one of three spins to the z direction for A > 1. 
The ground state is slightly different from the 120° state 
in the Ising case A > 1, but this deviation does not af- 
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Fig. 1. (Color online) (a) Dynamical spin correlation G(t) as a 
function of the Monte Carlo step t for the model (1) in the XY 
limit A = 0. (b) Spin relaxation time r as a function of tem- 
perature T at A = 0. The curve shows a fit by the BKT scaling 
t = aexp[6/(T — Tbkt) 1//2 ] ■ Tbkt an( i T are the estimates 
of the BKT transition temperature and the onset of the BKT 
critical region, respectively. See the text for details. The inset 
demonstrates the scaling plot for G(t). Below T, all the data are 
scaled well to a single universal function (we discard the data for 
t < 100). 



feet the long-time behavior of relaxation process. We set 
J = 1 and the Boltzmann constant fee = 1. 

In Fig. 1, we demonstrate how the NER method works 
in the XY limit (A = 0), as an example. We calculate 
the dynamical spin correlation function G(t) defined as 

1 



G(t) = 



(Si(t)-Si{0)), 



(2) 



where Si(t) denotes the spin configuration at site i and 
MC step t. As shown in Fig. 1(a), G{t) changes from an 
exponential decay in the high-T paramagnetic phase to 
a power-law decay in the low-T BKT phase. The BKT 
transition temperature Tbkt is determined by the di- 
vergence of the relaxation time r estimated from the 
high-T exponential behavior, G(t) ~ exp(— t/r). We em- 
ploy the scaling analysis by using G(t) = t~ p / (t / V) , 20 
which enables us to estimate r up to ~ 10 5 . The re- 
sults are plotted in Fig. 1(b). The divergent behav- 
ior of r at low T is well fitted by the BKT scaling 
r = g(T) = aexp[6/(T-T BK T) 1/2 ]- Here we choose the T 
range for the fit by monitoring the weighted residual de- 
fined by R w = ^ Y^ilin - 9(Ti)}/g(Ti)]\ where N T 
is the number of the data {T}: We fit the range of low- 
T data which gives i? w less than 0.002. The fit gives an 
estimate of the transition temperature Tbkt = 0.409(1), 
which is consistent with the previous estimate. 16 At the 
same time, the fitting procedure defines T = 0.440(10), 
below which r follows the BKT scaling. The region 
Tbkt < T < T represents the BKT critical region in 
which G(t) obeys the universal behavior [see the inset of 
Fig. 1(b)]. We confirm a similar universal scaling for all 
the following data. 

We study the critical behavior of r in this way for 
various values of the anisotropy A. The results for A < 1 
are shown in Fig. 2(a). As approaching the Heisenberg 
case with A —> 1, Tbkt decreases monotonically, while 
T first decreases but increases for A > 0.9: As a result, 
the width of the BKT critical region becomes wider as 
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Fig. 2. (Color online) (a) Critical behaviors of r for various values 
of the anisotropy A. (b) BKT scaling plot for the normalized re- 
laxation time f . The lines show the BKT scaling fit. (c) Relative 
width of the BKT critical region as a function of the anisotropy. 
(d) T dependence of r at the Heisenberg point A = 1. The curve 
shows the BKT fit, and the dashed line represents the fitting 
on the basis of the renormalization-group (RG) analysis for the 
non-linear cr model, r = C t [(T/B) x exp(B/T)] z . 21 



A — y 1. This is more clearly observed in Fig. 2(b), which 
plots the normalized relaxation time f = rja on the 
basis of the BKT scaling r = acxp[6/(T - T BK t) 1/2 ]- 
In fact, as shown in Fig. 2(c), the relative width of the 
critical region, ATbkt = (T — Tbkt)/Tbkt, appears to 
diverge logarithmically or more strongly with decreasing 
the anisotropy. 

As anticipated from the diverging ATbkt, t for the 
isotropic Heisenberg case can be fitted by the BKT scal- 
ing in the entire range of T calculated, as shown in 
Fig. 2(d). The fitting naively suggests that t diverges at 
T* = 0.282(4). Similar behavior was seen in the spin cor- 
relation length, for which it was argued that a crossover 
takes place from the BKT behavior to another T depen- 
dence and the correlation length does not diverge ex- 
cept for T = 0. 22 ' 23 Our data are consistent with such 
analyses as shown in Fig. 2(d). We return to this point 
later. Besides the crossover, the crucial observation here 
is that the isotropic case A = 1 looks quite singular since 
the apparent BKT critical behavior is observed in the 
divergently wide range of T. 

Let us further discuss the singular behavior as A — >• 1 
from the viewpoint of the vector spin chirality. We study 
the dynamical chiral correlation function defined as 

(3) 
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Fig. 3. (Color online) Relaxation of the vector chiral order pa- 
rameter (a) in the XY limit (A = 0.0) and (b) very close to 
the Heisenberg point (A = 0.9999). (c) Scaling plot of the relax- 
ation time of the chirality, t k . Solid lines for A < 0.99 represent 
the power-law fit r K oc (T — T c )~ zv ■ The dashed curve for the 
Heisenberg case A = 1 is the fit to the BKT criticality, i.e., 
t k oc exp[6/(T — T*) 1 / 2 ]. (d) Relaxation in the Heisenberg case. 



where the summation Ri runs over all the unit trian- 
gles, and k z is the z component of the vector chirality 
k z = ( Si x S 2 + 5 2 x S 3 + S 3 x Si ) z for each triangle of 
three spins Si, S2, and S3. The vector chirality exhibits a 
true long-range order in the anisotropic cases, and there- 
fore, wc expect that n{t) decays exponentially above a 
chiral transition temperature T c , while it approaches a 
nonzero constant below T c : At T = T c , n(t) shows a 
power-law decay. This is indeed the case as demonstrated 
in Figs. 3(a) and 3(b) for A = 0.0 and 0.9999, respec- 
tively. The results give the estimates T c = 0.413(1) for 
A = 0.0 16 and T c = 0.290(5) for A = 0.9999. It is note- 
worthy that the chiral transition is clearly discernible 
even for very small anisotropy of 0.01% (A = 0.9999). 

Similar to the analysis of r in Fig. 2, we examine the 
behavior of the relaxation time of the chirality, t k , by 
varying the anisotropy A. The results for A < 1 are sum- 
marized in Fig. 3(c). For A < 1, T K shows a power-law 
divergence t k oc (T — T C )~ Z1J ', with the same exponent 
zv ~ 1.9. At A = 0.999 the data show crossover from a 
BKT-likc behavior at high T to the power-law divergence 
near T c . In the isotropic case A = 1, r K shows a stronger 
divergence than the power law and is well fitted by the 
BKT scaling r K ex exp[b/{T-T*)^ 2 } with T* = 0.284(3). 
The relaxation process at A = 1 is shown in Fig. 3(d); 
k(£) exhibits a power-law decay even at much lower T 
than T* . The results seemingly suggest a BKT-type tran- 
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Fig. 4. (Color online) (a) Phase diagram for the anisotropic 
Heisenberg model (1) determined by the NER method. For the 
XY anisotropy (A < 1), the chiral and BKT transition temper- 
atures are shown by diamonds and crosses, respectively. For the 
Ising anisotropy (A > 1), two BKT transition temperatures as to 
S 2 and (S x , S y ) components are plotted by triangles and circles, 
respectively. The lines are guides for the eye. (b) Chiral and BKT 
transition temperatures as a function of 1/| ln(l — A)| for A < 1. 
The data are fitted by 1/| ln(l — A) I". For comparison, a recent 
estimate of Z2 vortex transition temperature T v is shown. 25 



sition in the vector chiral degree of freedom at T* , as seen 
at T* in the spin sector [Fig. 2(d)]. These behaviors in 
the vector chirality also illuminate a singularity of the 
isotropic Heisenberg case. We will discuss T* and T* in 
comparison with the Z2 vortex transition temperature 
later. 

Collecting the data of the relaxation of spin and vec- 
tor chirality, we map out the finite-T phase diagram as 
a function of the anisotropy A. We also study the Ising 
anisotropic cases (A > 1), where two different BKT tran- 
sitions take place; one is for the z component of spin S z , 
and the other is for the xy components S x , S v . 15 The 
latter transition accompanies a quasi long-range order- 
ing of the vector chirality, and occurs at a lower T than 
the former. 18 The two BKT transition temperatures are 
determined from the dynamical spin correlation function 
for the corresponding spin component [cf. Eq. (2)]. 

Figure 4(a) summarizes our phase diagram around the 
isotropic Heisenberg point A ~ 1. In the XY anisotropic 
region A < 1, the chiral transition always occurs at a 
slightly higher T than the BKT transition. 16 Both two 
transition temperatures decrease more rapidly as A — > 1. 
Similar decrease is observed in the two BKT transitions 
when approaching from the Ising anisotropic case A > 1. 
Especially, the BKT transition temperature for the xy 
components shows a nonmonotonic A dependence. Con- 
sequently, the phase diagram exhibits a sharp "V shape" , 
which illuminates the singularity of the Heisenberg case 
A = 1. To our knowledge, this peculiar form in the very 
vicinity of A = 1 has not been elucidated before. 24 

The question is the fate of the transition temperatures 
as A — > 1, in particular, their relation to the Z2 vortex 
transition predicted for A = 1. Figure 4(b) shows the 
asymptotic behavior of the chiral and BKT transition 
temperatures in comparison with a recent estimate of Z2 
vortex transition temperature T v . 25 We plot the data as 
a function of 1/| In (1 — A)|, with considering an analyti- 
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cal argument for the square lattice model which predicts 
?bkt oc 1/| In (1 — A) | on the basis of independent vortex 
pair picture. 26 Surprisingly, both T c and Tbkt decrease 
faster than oc 1/| ln(l — A)|, suggesting that they will be 
well below T v and finally approach zero as A — > 1. (The 
fitting in the figure shows 1/| ln(l — A)| Q with a < 1 as 
a guide.) From these observations, we conclude that the 
phase boundaries of the chiral and BKT transitions show 
the highly-singular "V-shape" around the isotropic point 
A = 1, and the Z 2 vortex transition is isolated from the 
conventional phase transitions. 

We note that the Z 2 vortex transition temperature 
T v = 0.285(5) 25 coincides with the apparent 'transition 
temperatures' T* = 0.282(4) estimated from the BKT 
fitting of t [Fig. 2(d)] and T* = 0.284(3) similarly ob- 
tained for t k [Fig. 3(c)]. Since the relaxation time r and 
t k exceed 10 6 , it is hard to trace a crossover from the 
BKT scaling to another behavior, if any, within the acces- 
sible system size [see Fig. 2(d)]. Although dynamics of Z 2 
vortices may cause these diverging behaviors, their re- 
lation is not clear. Nonetheless, as discussed above, since 
we expect a finite spin correlation length for T > at 
A = 1, it is natural to consider Tbkt — ► as A — >• 1, not 
Tbkt — ^ T*. The situation is not so clear for the vector 
chirality, but it is plausible that T c also goes to zero as 
A —> 1, because T c coincides with Tbkt for S x , S v in the 
Ising case A > 1, which should go to zero. The asymptotic 
behaviors in Fig. 4(b) support our consideration. 

Finally, let us discuss the relevance of our results to 
experiments. The peculiar "V shape" phase diagram in 
Fig. 4 indicates that the anisotropy is a relevant pertur- 
bation to the isotropic Hcisenberg point A = 1 in the 
sense not only that it triggers the finite-T transitions 
but also that the induced transition temperatures grow 
in a very singular fashion against the anisotropy. There- 
fore, in the triangular antiferromagnets, the isotropic 
Heisenberg case is special, rather isolated; the pristine 
property of the isotropic point including the Z 2 vor- 
tex transition is hardly accessible in real materials in 
which anisotropy exists inevitably. For example, NaCrC>2 
is known to have a small but finite anisotropy in the g 
value about 0.25%, 27 and NiGa2S4 has a rather substan- 
tial anisotropy of ~ 3%. 14 Although there might be a dif- 
ference between the exchange anisotropy and the single- 
ion anisotropy, 28 our results strongly suggest that in real 
materials the conventional BKT transition dominates the 
critical behavior of two-dimensional spin fluctuations, in- 
stead of the unconventional Z 2 vortex transition. The 
anomalous enhancement of the spin relaxation time will 
be understood by the divergently-enlarged BKT critical 
region near the isotropic case. 

A potential "smoking gun" for the relevance of 
anisotropy is critical behavior of the relaxation time of 
the vector chirality. As shown in Fig. 3(c), the power-law 
criticality is observed in the presence of the anisotropy, 
whereas the BKT-type divergence dominates at the 
isotropic point. It is difficult to probe the chirality experi- 
mentally, but we note that recently the chiral phase tran- 
sition was detected by polarized neutron scattering. 29 
The development along this direction is highly desired. 

In summary, by using the nonequilibrium relaxation 



Monte Carlo method, we have studied the relaxation 
time of spin and vector chirality in the anisotropic clas- 
sical Hcisenberg model on the triangular lattice. For the 
spin relaxation time, we have revealed that the BKT crit- 
ical region becomes divergently wide as the anisotropy 
decreases. For the vector chirality, the relaxation time 
exhibits a power-law divergence, whereas at the isotropic 
Hcisenberg point, it shows an apparent BKT criticality in 
a wide range of temperatures. We have also obtained the 
precise phase diagram with its singular "V shape" around 
the isotropic Heisenberg point, which uncovers the sin- 
gular nature of the isotropic point. Our findings will be 
important for understanding of the puzzling experimen- 
tal results in the triangular antiferromagnets as well as 
of unsettled theoretical issue on the relation between the 
conventional transitions and the Z 2 vortex transition. 
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